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Abstract 

Randomly-forced fluid flow in the presence of scale-unselective dissipation de- 
velops mean currents following topographic contours. Known mechanisms based 
on the scale-selective action of damping processes are not at work in this situa- 
tion. Coarse-graining reveals that the phenomenon is a kind of noise-rectification 
mechanism, in which lack of detailed balance and the symmetry-breaking pro- 
vided by topography play an important role. 
PACS numbers: 05.40.-a, 47.90. -Fa, 92. 90. -Fx 
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In the last decades much effort has been devoted towards understanding ffuid motion in sit- 
uations in which vertical velocities are small and slaved to the horizontal motion. Under these 
circumstances the ffow can be described in terms of two horizontal coordinates, the vertical 
depth of the fluid becoming a dependent variable. The fluid displays many of the unique prop- 
erties of two-dimensional turbulence, but some of the aspects of three-dimensional dynamics 
are still essential. Among the several quasi two-dimensional dynamics considered, barotropic 
quasigeostrophic turbulence has focussed most of the interest. The reason for such interest 
lies on the relevance of this dynamics as a model to understand planetary atmospheres and 
ocean general circulations as well as to describe some plasma physics phenomena In the 

first case the two dimensionality of the flow is induced by the Coriolis force, whereas in the 
second it may arise from the action of magnetic fields. In the geophysical context, the topog- 
raphy over which the layer of fluid flows is the ingredient introducing a difference with respect 
to purely two dimensional turbulence. 

The first works on topographic turbulence focused on its statistical properties in the ab- 
sence of dissipation and forcing 0. These studies estabhshed a tendency of the flow to reach 
a maximum entropy (Gibbs) state characterized by the existence of stationary mean currents 
following isobaths. Mean currents associated with topography are however not restricted to this 
inviscid and unforced shown by |^ , where viscosity was included. Statistical- mechanics 

equilibrium arguments can not be invoked in this situation of decaying turbulence to explain the 
appearance of mean currents correlated to topographic features. The explanation put forward 
was that viscous damping, due to its stronger action at the smallest scales, dissipates enstro- 
phy much faster than energy, a quantity more concentrated at large scales in twodimensional 
turbulence settings. As a consequence of this scale-selective behavior, the decaying turbulent 
state would be the one with the smaller enstrophy compatible with the (slowly decreasing in 
time) instantaneous value of the energy. It turns out that these minimum enstrophy states are 
closely related to the equilibrium maximum entropy ones and, as them, present distinct mean 
currents following isobaths. The natural tendency of topographic turbulence to generate cur- 
rents along isobaths, under different situations, aroused a significant interest in the geophysical 
community. Specifically, studies were addressed to investigate the role played by these currents 



on the general circulation of the world's ocean pHTT 
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It has been recently shown numerically [T^ that mean currents following isobaths appear 
also in situations in which a scale-unselective dissipation, Rayleigh friction modeling an Eckman 
bottom drag, is used. This fact was observed in numerical simulations in which random forcing 
(also of a scale unselective nature) was present. Thus the state obtained is not one of decaying 
turbulence but a nonequilibrium statistically steady state. Statistical properties are not far 
from those of a generalized canonical equilibrium, but this is just a numerical coincidence, since 
neither the maximum entropy nor the minimum enstrophy mechanisms, the ones invoked in 
the previous situations, are here at work. Instead ||12| interpreted these currents as a noise 
rectification phenomenon ||13[: the mean currents were generated by nonlinearity and sustained 
by noise, with topography providing the symmetry-breaking ingredient giving a particular sense 
to the flow. Later it was theoretically shown and numerically confirmed that noise rectifica- 
tion also appears in randomly forced topographic turbulence when a scale selective damping 
(viscosity) is considered 0,|l5l. However, the detailed understanding of the noise- rectification 
mechanism when a scale-unselective damping is present is still an open question. Addressing 
this question as well as to provide a rigorous theoretical framework to explain the numerical 
results shown in are the main motivation of this Letter. 

In the quasigeostrophic approximation, the motion of a single layer of fluid in a rotat- 
ing frame or planet can be described in terms of the horizontal components of the velocity 
{u{x.),v{x.)) that can also be written in terms of a streamfunction verifying: 

u = -^ v = ^ (1) 

dy ' dx 

The streamfunction '?/'(x, t), with x = {x,y), is governed by the dynamics @: 

— + X ^,V^ij + h = -eVV + F . (2) 
at L J 

e is the friction parameter, modelling bottom friction, -F(x, t) is any kind of relative- vorticity 
external forcing, and h = fAH/Ho, with / the Coriolis parameter, Hq the mean depth, and 
AH{x) the local deviation from the mean depth. A is a bookkeeping parameter introduced to 
allow perturbative expansions in the interaction term. The physical case corresponds to A = 1. 
The Poisson bracket or Jacobian is defined as 

\A,B] = ——-——. (3) 
dx dy dx dy 
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Equation @ gives the time evolution of relative vorticity subjected to forcing and dissipa- 
tion. The form used by the dissipation term acts with the same strength at all spatial scales, 
in contrast with viscosity terms of the form V^ip, excluded from our model (0), which would 
dissipate better the smaller scales. A convenient choice of F, able to model a variety of pro- 
cesses, is to assume it to be a Gaussian stochastic process with zero mean, and with a Fourier 
transform F]^{u) having the following two-point correlation function: 

Fk(cu)Fk'(^')) = Dk-y6{k + k')6{uj + uj') . (4) 

Here, k = (fc^, ky), and k = |k|. Thus the process is white in time but has power-law correlations 
in space. Stochastic forcing has been used in fluid dynamics problems to model stirring forces 
||T6|, short scale instabilities thermal noise or processes below the resolution of 



computer models among others pO 



Guided by the results in [Q, we attribute the average currents to a large-scale rectification 
of the small-scale fluctuations introduced by the noise term. In consequence we resort to a 
coarse-graining procedure to investigate how the dynamics of long-wavelength modes in is 
affected by the small scales. For our problem it is convenient to use the Fourier components 
of the streamfunction ■j/'koj or equivalently the relative vorticity = —k'^i'kuj- This variable 
satisfies: 

Ckoj = Ckw-^ki^ ~^ 

p,q,n,n' 

where the interaction coefficient is: 

^kpq = (PxQy - Pyqx)P''^ Kp+ci , (6) 

the bare propagator is: 

GL = H^ + e)-^ , (7) 

and the sum is restricted by k = p + q and uj = Q + Q'. p = {px,Py), P = |p|, and similar 
expressions hold for q. The wavenumbers are restricted to < k < ko, with ko an upper 



cut-off. Following the method in Ref. one can eliminate the modes CkL "with k in the 



shell k^e^^ < k < ko and substitute their expressions into the equations for the remaining 
low-wavenumber modes with < A; < k^e"^ . The parameter 5 measures thus the width of 
the band of eliminated wavenumbers. To second order in A, the resulting equation of motion 
for the modes CkL, written in terms of the large-scale streamfunction ip^{'x,t) is: 



dt 

- eVV^ - 9^^h< + Z/'VV^ + F' , (8) 



where 



16(27r)2e2' ^""^ 

\^DS2{y + 2)5 

9(>^^D,6,e,y) = ^^^^^^ . (10) 

F'(x, t) is an effective noise which turns out to be also a Gaussian process but with mean 
value and correlations given by: 

^ ' ' 32(27r)2e2 ' ^ ' 



Dk-y6{k + k')6{uj + uj') (12) 

5*2 is the length of the unit circle: 2tt. Equations (p|)-(p!^ are the main result in this Letter. 
They give the dynamics of long wavelength modes ip^ with < k < k^e'^ . They are valid for 
small A or, when A ~ 1, for small width 5 of the elimination band. 

The elimination of the small scales leads, as expected from physical grounds, to the appear- 
ance of an effective viscosity v' at large scales. Depending on the sign of ?/, this viscous term 
can be destabilizing. One should keep in mind however that Eq.(^ is only valid at large scales, 
so that such small- wavelength instabilities, when formally present, would be avoided by proper 
choice of the wavenumber cut-off. More importantly, a new term depending on the topography 
gV^h^ has been generated by the small scales. Another similar term is contained in the mean 
value of the effective large-scale noise F'. Both terms have the effect of pushing the large-scale 



5 



motion towards a state of flow following the isolevels of bottom perturbations h'^. The energy in 
this preferred state is determined by the function g{X, D, 6, e, y) which measures the influence of 
the different terms of the dynamics (nonlinearity, noise, friction). From relation ( ]T0| ) it is clear 
that while nonlinearities and noise increase the intensity of the mean currents in the preferred 
state, high values of the friction parameter would reduce them. 

As in the situation with scale selective damping [ p^p3| , (|1^) and (|Tl]) exhibit a characteristic 
dependence on the exponent y of the random-forcing spectral power-law. In the present case, 
the dependence with y + 2 implies that the directed currents reverse sign as y crosses the 
value —2, and that they vanish ii y = —2. It is easy to check that the conditions for detailed 



balance between the fluctuation input and dissipation p2| are satisfied precisely for y = —2 
(in general, the condition is y = —n, where n is the exponent of V in the dissipation term). 
In this situation the steady-state probability distribution for can be found exactly and is 
independent of topography. Thus the generation of flow along isobaths is a consequence of lack 



of detailed balance, a general feature of noise-rectifying mechanisms |13 



Concluding, relations (|8|)-([T^ show that the origin of average circulation patterns in quasi- 
geostrophic turbulence over topography is related to nonhnearity and lack of detailed balance. 
The present approach provides a theoretical basis for the numerical observations presented in 
|T^, for which mechanisms based on the scale-selective action of damping can not be applied. 
Nonlinear terms couple the dynamics of small scales to the large ones and provide the mecha- 
nism for energy transfer from the fluctuating component of the spectrum to the mean one. This 



mean spectral component, absent in purely two-dimensional turbulence ||23|, is controlled by 
the shape of the bottom boundary, which breaks the isotropy of the system, and characterizes 
the structure of the flow pattern. 
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